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Abstract 

A study of the thermal properties of two-dimensional topological lat- 
tice models is presented. This work is relevant to assess the usefulness of 
these systems as a quantum memory. For our purposes, we use the topo- 
logical mutual information /topo as a "topological order parameter". For 
Abelian models, we show how /topo depends on the thermal topological 
charge probability distribution. More generally, we present a conjecture 
that /topo can (asymptotically) be written as a KuUback-Leitner distance 
between this probability distribution and that induced by the quantum 
dimensions of the model at hand. We also explain why /topo is more suit- 
able for our purposes than the more familiar entanglement entropy Stopo. 
A scaling law, encoding the interplay of volume and temperature effects, 
as well as different limit procedures, are derived in detail. A non-Abelian 
model is next analysed and similar results are found. Finally, we also 
consider, in the case of a one-plaquette toric code, an environment model 
giving rise to a simulation of thermal effects in time. 

1 Introduction 

Quantum mechanics has significantly shaped our current understanding of con- 
densed matter systems. It has proven to be deeply insightful when aiming to 
explain, say, transport of charged (quasi-) particles in semiconductors, mag- 
netism in metallic alloys or cohesive properties of solids [T . But when it comes 
to macroscopic systems that exhibit a highly non-classical behavior, many fun- 
damental issues are still poorly understood. An example of such systems is a 
fractional quantum Hall effect sample (FQHE), i.e. a specific two-dimensional 
electron gas subject to a strong perpendicular magnetic field. There the trans- 



verse conductivity appears in plateaux at fractional values of the filling factor, 
much in contrast to what one would expect classically Besides their fun- 
damental interest, the exotic phases exhibited by these systems are important 
in that they may lead to new technological applications. Indeed, they could 
allow for intrinsically fault-tolerant quantum computation [31 |J . These phases 
are not separated by a symmetry whose absence or presence can be detected by 
a local order parameter. Rather, it has been realized that they are associated 
with topological order. 

In the effort towards understanding this notion, an important contribution 
has been the research towards lattice spin models. Certain such systems are 
exactly solvable [21 HI El E] , in the sense that their low energy sectors can be 
analytically determined. Although the corresponding interactions are local and 
frustration- free, they possess the essential properties of topologically ordered 
systems: (i) when the system is defined on a surface with non-trivial topology, 
the ground state has a non-local degeneracy that cannot be detected locally, 
(ii) excitations have exotic statistics. These features turn these models into 
an interesting alternative to FQHE systems, when it comes to study topologi- 
cal order. Initially, topological spin lattice systems were proposed as a reliable 
quantum memory [3 dill] and some non-Abelian versions were shown to allow 
for universal fault-tolerant quantum computation [101. Quantum information 
would be encoded in non-local degrees of freedom of such systems and would 
therefore be immune to local perturbations. Fault-tolerant quantum computa- 
tion could be performed by creating excitations (initialization), braiding them 
(unitary evolution), and fusing them back together (read-out) [^ 

In this work, which is partly an extension of |12j . we focus on an entropic 
order parameter, the topological mutual information, /topo, and use it to study 
how topological matter behaves in the presence of temperature. For concrete- 
ness, we restrict our analysis to two paradigmatic models: the toric code and 
the 0(33) "superconductor" [31|T3]. We will see that at any fixed finite tem- 
perature, /topo is non-zero only when the size of the system is finite. Next, for 
fixed finite system size, there is always a temperature regime where the order 
parameter assumes a constant value before dropping to zero. A similar behav- 
ior is observed when the size of the system is increased for a fixed value of the 
temperature. Importantly, we exhibit a scaling relation that tells how much the 
temperature should be decreased to compensate for an increase of the size of the 
system, if the system is to remain topologically ordered. The exact behaviour 
of topological ordered systems with respect to temperature has received much 
attention recently. From a static point of view, the first results indicating the 
fragility of the 2D toric code against temperature can be found in [TH [15] . More 
recently, the authors of [HIITT] have shown the fragility of this system during the 
thermalization dynamics. In |181 119] a more general approach is taken to show 
the impossibility of having a self correcting quantum memory if one restricts 
the search to 2D stabilizer codes. Apart from providing exact scaling relations, 
our work is the first to deal with the non-Abelian situation. A different new ap- 
proach in this respect can be found in ^20j. Though /topo is initially introduced 
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as the natural generalization of the topological entropy S'topo of Kitaev-Preskill 
[2T] and Levin- Wen [22' to the thermal situation, the usefulness of this quantity 
will become clear throughout the paper. In particular, at least for the toric 
code, we explicitly relate /topo with the probability distribution of topological 
sectors in a region, connecting in this way the value of /topo with the ability of a 
system to perform quantum computing tasks at finite temperature. Moreover, 
we show how /topo allows to distinguish between quantum double models which 
share the same value of S'topo at zero temperature. 

The paper is structured as follows. In Section [2] we introduce /topo and 
discuss some of its properties. In Section |3] we compute /topo for the toric code 
and show explicitly its dependence on size and temperature. We also explain 
how to simulate the behavior of /topo in a one-paquette toric code with current 
technology. In Section |4] the behavior of /topo for the non-Abelian case is 
studied. The general formalism is presented and the special case of 0(83) is 
given in detail. We provide numerical evidence of the fact that /topo depends 
on size and temperature exactly as for the toric code. Finally, in Section |5] we 
discuss the conclusions and implications of our work. 

2 Topological mutual information 

Our analysis is based on the idea that constant corrections to area laws are 
typical signatures of topological order. Let us have a closer look at this property. 
Consider a bipartition /? : /?° of a given system in a pure state of its ground 
subspace, and assume that the von Neumann entropy of the reduced density 
operator of /?, = —tTpRln pj^ satisfies: 



where a' is a constant, \dR\ denotes the size of the boundary of R and where e 
tends to zero when \dR\ tends to infinity. As discussed in [51] , systems 

with a non-zero value of the topological entropy, 7', are topologically ordered. 
Indeed, 7' is related to the total quantum dimension V of the anyonic model 
describing the excitations: 



where dq is the quantum dimension associated with anyon type q |22l 121] . The 
trivial case with T) = 1, i.e. 7' = 0, corresponds to non-topological models, 
where the only contribution to V comes from the vacuum. 

Let us consider a system defined on a closed surface E. In the remainder of 
this paper, E will either be a torus or a sphere. It was shown in [55] |5T| that 
7' can be expressed as a linear combination of entropies of regions of E. For 
example, if S is divided into four regions, as indicated in Fig. [T| then ^21] 

— Sa + Sb + Sc — Sab — Sac — Sbc + Sabc- (3) 



SR^a'\dR\-y + e{\dR\), 



(1) 




(2) 
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Figure 1: Division of a torus or a sphere into four regions. 



Actually, there is much freedom in constructing linear combinations of entropies 
that are topologically invariant, i.e. invariant under local deformations of the 
boundaries. It is reasonable to require that the regions A, B, C be treated on 
the same footing, and thus write such a linear combination as /topo = ci^iSA + 

SB + Sc) + a^{SAB + SAC + SBc) + a^SABC + b^iSAD + SBD + ScD) + b^{SABD + 

SACD+SBCD)+b^ Sabcd+xSd- Consider a deformation of the frontier between 
regions C and D away from any triple point. Since the deformation is well inside 
the region CD, we have AScd = ^Sacd = ^Sbcd = ^Sabcd = 0. Since it 
is local between the regions C and D, AS'c — /S.Sac = ^Sbc = ^Sabc and 
A5'r, = AS'Ai, = ASbd = ASabd- Therefore 

A/t„p„ = (fli + + a^)ASc + {2b^ + b^ + x)ASd. (4) 

So, in order to get a topological invariant, we must have that a^+2a^+a^ — 2b^ + 
+ X = 0. Now let us consider a triple point deformation, at the intersection 
between the regions B, C and D, say. Reasoning as before, we get 

A/topo = (a' + a^){ASB + ASc) + (a^ + a')A5Bc 



+ib^ + b^){ASBD + AScd) + {b" + x)ASd. (5) 

We thus get stronger conditions: — —a? = and b^ — —b^ = —x. The 
entropy of the total system Sabcd is irrelevant, as expected; it is invariant 
under boundary deformations. Hence, there is no constraint on b^. For = 
l,b^ = b'^ = X — 0, we recover the topological entropy, 7', defined in [21]. In 
the following, we work with the choice = —b^ = —b^ = x = 1. This choice 
yields 

hopo ^ Ia + Ib + Ic - Iab - I AC - Ibc + Iabc, (6) 

which amounts to replace the von Neumann entropies appearing in the definition 
of the topological entropy [21] by quantum mutual information. is defined 
as Sr + Sr^ — Srur^-) 

At finite temperature, the von Neumann entropy of a region i?, Sr, is not a 
measure of correlations between R and the rest of the system, as it is in the pure 
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state case. In contrast, the quantum mutual information, Ifi, still is. Moreover, 
Sn does not obey an area law of the form (jlj anymore, whereas for the lattice 
systems we are going to study, the mutual information still does That is, 
the properties of the von Neumann entropy that make 7' a topological order 
parameter at zero temperature are no longer valid at finite temperature. But 
they still hold for the quantum mutual information. This is why we choose to 
work here with /topo = 7 instead of 7'. In Section 3.5 we will further discuss 
the behaviour of 7'. 

We close this section by showing that any linear combination of entropies 
that is topologically invariant should, in general, involve a division of the surface 
E into at least four regions. It will be enough to consider the case where 
the whole system is in a pure state. We would like to construct a particular 
linear combination of entropies that isolates 7 from the area part. It is clear 
that partitioning the surface into a region R and its complement R'^ cannot 
provide such a quantity. Indeed, in this case Sr — Sr^, Srur^ — and dR 
and I7I have a common fate. Three regions A, B and C are not sufficient 
neither. This is easily shown considering a system in a pure state. Let Iat^b 
and Iab denote respectively the length of the boundary between region A and 
region C, region B and region C, and region A and region B, and let nA 
denote the number of connected pieces that make region A, ub and uab are 
defined likewise. From the relations Sa + Sb = + Ib + Iab) — (nA + nB)j, 
Sa-Sb ^ a{lA ~Ib) + {riA - ?^B)7 and Sab = Sc = a{lA + h) ~ nABl, we see 
that it is impossible to construct a linear combination of entropies that cancels 
all boundary contributions and leaves only the topological contributions. Thus, 
four is the minimal number of pieces required in order to partition S in such a 
way that it gives a boundary-independent quantity in a non-trivial way. 



3 The toric code 

The toric code is a simple topological model, with a Hamiltonian that can be 
diagonalized exactly [5]. It is called a 'code' because it is a quantum error 
correcting code; two logical qubits are encoded in the physical system. An 
optical implementation of a four-qubit toric code is possible since its ground 
state is a GHZ state [25], and proposals to create the ground state of large 
toric codes, as well as elementary excitations in optical lattices are described in 
[25] . In this section, we study how /topo behaves as the (inverse) temperature, 
/3, and the size of the system, L, are varied. A similar calculation has been 
presented in [TS]. We nevertheless present our alternative approach in details 
because it differs in two important respects. First, it is 7' which is used in 
[15] as a topological order parameter. Second, the simphcity of our alternative 
analysis has allowed us to get analytically more general as well as new results 
such as the scaling laws discussed in Section |3.3[ Moreover, this calculation 

^Note that under quite general assumptions, a weak form of area law always holds at finite 
temperature: Iji < a\dR\, for some constant a |24| . 



5 



helps understand better how to compute the topological mutual information for 
non-Abelian models. 



3.1 Spectrum 

Let us consider a torus tiled into L x L square plaquettes and associate a two- 
level system (qubit) with each edge of the obtained lattice. We assume that 
these qubits interact through the Hamiltonian 

H=-JY,Bp-J'J2As, (7) 

p s 

where the index p (resp. s) runs over all plaquettes (resp. vertices) of the tiling. 
The operator Bp involves all the spins surrounding the plaquette p, while the 
operator As involves all the spins with one end at s. They are defined as follows: 



A, 



(8) 



and represented on Fig|2] The coupling constants J and J' will be chosen to 
be both equal to 1, for the sake of simplicity. But all our analysis can be 
straightforwardly generalised to arbitrary values of J and J'. 



Figure 2: (Left) Pictural representation of the Bp and the As operators. (Right) 
An eigenstate of the toric code lattice. Elementary excitations can be either of 
flux type, depicted by squares in plaquettes, or of charge type depicted by hollow 
dots on vertices. The strings along the torus correspond to Wilson loops. 

Due to the topology of the torus, we have that 

IlpBp = UsAs = 1. (9) 

This constraint means that the number of excited plaquettes (resp. excited ver- 
tices) is always even. The Hamiltonian _ff is a sum of local terms all commuting 
with each other. We observe that although this property facilitates the diago- 
nalization of H, it is not sufficient to guarantee that it is easy to solve. What 
makes H exactly diagonalizable is that, as it turns out, its eigenstates of H 
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can have arbitrary eigenvalues of the operators Bp and As , up to the constraint 
([9|. Therefore, each eigenstate of H is given by a triple of "quantum numbers": 
|(/), c, w), see Fig. [2j A pattern </) denotes the position of all plaquette or "flux- 
type" excitations. Another pattern, c, indicates the position of all vertex or 
"charge-type" excitations. Finally, w indexes the degeneracy of the state for a 
fixed configuration of defects. This quantum number is made of two bits, wi 
and W2, that label the values of the integrals of motion of z-operators around 
non-contractible loops on the torus [3] (Wilson loops). We have 

H\^, c, w) - {Eo + m + 2\c\) 10, c, w), (10) 

where Eq = —2L^ is the ground state energy and \4>\ (resp. |c|) denotes the 
number of flux excitations (resp. charge excitations) of the pattern 4> (resp. c). 
The eigenvalues of H satisfy £"„ — £'„+i = —4 and range between —2L^ and 
2L^. Also, if Pi denotes the projector onto the sector of energy Ei — Eq + 4:i 
and di = tr Pi denotes its dimension, we have that 

One can check that di = 2^^ . 
3.2 von Neumann entropy 

In the following, we shall consider a situation in which the system is immersed in 

a bath at inverse temperature (3 and is let to thermalize. Since H = X]i=o ^i^ii 
the partition function of this model reads 

Z{P, L) = tr e-P" = ^ e-^^^d^, (12) 

i=0 

This series can be easily summed up (see [H] or details in Appendix |A|. We get 

Z(/3,L) = ((2 cosh /3)^' -I- (2sinh/3)^')2. (13) 

The thermal state of Eq.([7]) reads 

p,^^e-P"/Z{(3,L). (14) 

The von Neumann entropy, Stot — — trpth In Pth, of the whole torus is then easily 
derived from the partition function thanks to the identity 

^to. = -^|mL) + lnZ(/3,L). (15) 

We now compute the von Neumann entropy of a connected region i? C S, 
Sr. a couple of observations about the reduced state, pa, allows to get an 
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analytic expression for Sr. First, we consider a fixed eigenstate c, w) of the 
Hamiltonian Q. The reduced state pR{<j),c,w) = tr }ij(f>,c,w){(l),c,w\ does not 
depend on w if i? is contractible; homologically non-trivial loops are necessary 
to measure w. Also, two states {(p,c,w) ^ {(j)' ,c/ ,w') are orthogonal whenever 
{(I^R: Cr) ^ {4>'r^ c^) since they can be discriminated by measuring operators 
or Bp operators having support on R. Next, it is useful to distinguish three 
kinds of plaquette excitations: those with support fully on R, (j>R, those with 
support fully on R^ 4>r^, and the others, 4>dR- Similarly, we divide vertex 
excitations into three kinds: cr,cr^,cqr. In order to lighten the notations, we 
use the symbol q to label configurations of defects, both plaquette and vertex, 
i.e. q = {4>, c). Crucially, excitations of the dR type can be driven inside R^ by 
application of Pauli operators acting on links of Rc, that is: 

for some unitary operator U'^ and some configuration of defects in R,., q^ . 
This property is illustrated in Fig. [3j 




Figure 3: Solid edges belong to a region R, while dashed edges belong to R^ = 
Y\R. A plaquette excitation on the right dR plaquette can be driven inside Rc 
applying cr^ on the edge separating R from Rc- 

All excitations inside R can be fused into a single excitation q]^ located on a 
sit^ using a unitary whose support is fully in R and similarly for Rc. In turn, 
the state I?)?, 9)?^, can be created from a ground state \£,,w) by application 
of Pauli operators along strings connecting the site where q]^ is located to the 
site where q\ is located. So, 

|q, w) = C/flXqi?) ® Ur^ (qfl, , qofl) le, w) (17) 

for some unitary operators Ur{(\r) and Ur^ {^.Ra 7 qai?) • We are now in a position 
to characterise pr and compute Sr. The thermal state of the toric code can be 
written explicitly as 

g-/3(Bo+A£;|q|) 

pth = 2^ — ziR L) — (18) 

where AE = 2 is the energy associated with a single excitation (plaquette or 
vertex). Therefore, for a non-contractible region R, 

PR^Y. C{qR)trRAUR{ciR)\i,w){^,w\UR{^R)^l (19) 

^As in [3, we call a site a combination of a vertex and an adjacent plaquette. 
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where 



C(qfl)= e-''(^''+2i'i«i+2|q«J+2|qaHi)/z(/3,L). 



Note that 4C(qjj) is the marginal probability of a configuration of defects q^. 
The decomposition ^19\ allows to compute Sr. Indeed, from the identity 



5(0 A,p,) = - ^AanA,: + ^A,5(p,), (20) 

i i i 

we find that 

Sr = 5f - E ln(4C(q«)), (21) 

where 5^** is the von Neumann entropy of the region R when the system is in a 
pure ground state |C, w): S'g = {\dR\ - 1) ln2 P^. 

It turns out that the sums appearing in Eq.(21 1 can be carried out exactly 
(see details in Appendix |A| . The result is that the entropy of a region can be 
expressed as 

Sr = 5f + V{p, NpiR),L) + V{p, N,iR),L), (22) 

with 

V{P, N,L)^Nln{l + e-'P)+ ln(l + 0^') + ^-y^ + ^ 

2iiTW) + ^ ) 



(l_gA-)(l-ffL -AT) 

2{rTW) — ^"^^"^ ^^^^ 

where — tanh /?. 

As we can see, the entropy of a region separates neatly into a pure state 
contribution and a finite temperature contribution. The first only involves the 
area of the region, while the second depends on its volume {Np{R) and N^{R)). 
By increasing the temperature, we pass from an area law to a volume law, as ex- 
pected. The calculations are almost identical when R is only semi-contractible. 
The only difference is that pn will depend on only one of the Wilson loops, 
and the entanglement entropy picks a — In 2 additive correction. When R is 
completely contractible, the correction is — In 4. 
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3.3 Dependence of hopo on size and temperature 



Using Eq.(p3|, we have plotted the topological mutual information ^ for tori 
of various sizes, see Fig. |4] For a fixed torus size, we observe that there is a 
region of values of p such that /topo is non-zero and stationary. Then, as the 
temperature is increased, /topo smoothly vanishes. Interestingly, the transition 
does not become more abrupt when the size of the system is increased. Rather, 
the curves displayed are similar. 



topo 




P 
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Figure 4: Topological mutual information as a function of f3 for tori of size 
4^ X 4^, 4^ X 4^, 4^ X and 4^ x 4^. The size of the region A, as indicated in 
Fig. [l]is kxk^ (L/4) X (L/4). 



From Fig. |4] we also see that, when the system size increases, the temper- 
ature at which the transition occurs decreases. Actually, in the limit of large 
codes, the transition temperature vanishes. When systems described by a local 
hamiltonian are left in thermal equilibrium, their mutual information is bounded 
by a constant times their area |24j . As it turns out, the systems we are studying 
obey a strict area law, i. e. the quantum mutual information between a region 
R and the rest of the sytem satisfies: 



In^a{(3,L)\dR\-^{P,R,L). 
For the toric code, we can compute 
aoo(/3)- lim lR/\dR\, 

>oo 



(24) 



7oo(/3)= lim (iR^a^mdRD- 

\dR\~*ca 



Interestingly, we find slightly different results for 7oo(/3), depending on how 
this limit is taken. One possibility is to first consider the limit for L oo, and 
then let the size of the region R grow. One finds: 



7(/3,i?,oo) = 21n2 

H In h — 



QN,iR) i_0N,{R) 
— ; In 
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(2^) 

where 9 = tanh/3, and where Np{R) = — Np{Rc) denotes the number of 
plaquettes within a region R and at its border. N^,{R) is defined Ukewise. 



Eq. (25 1 lends itself to a simple interpretation. Let pe denote the probability 
that a plaquette or a site is excited. The mean energy of the system reads 
(H) ^ Eo + 4peL^. From {H)Z{(3,L) = -dZ{P, L)/df3, we find that in the 
limit L — > oo, = (1 — ^)/2. On another hand, ignoring total anyonic charge 
conservation, the probability that a region R containing Np{R) plaquettes has 
an even number of excited plaquettes reads 

(i?) ^I'^j^ (^'^.^^)iplH-Pey)il-Pef^^''^-' = ^(l+0^''(^)). (26) 
2=0 V * / 

One can similarly calculate p°/™{R), the probability that the region R con- 
tains an even number of excited vertices. One find the same expression with 
Np{R) replaced by N^,{R). So the asymptotic limit of the topological mutual 
information can be rewritten as 

7oo(/3) = 21n2 - /i2(pr"(^)) - h^ipT'^m, (27) 

where h2{x) — —a; In a; — (1 — a:)ln(l — x) is the Shannon entropy of a binary 
outcome probability distribution, or equally simply as 

7oo(/3) = -i?({Pfl(g)}||K(g)}), (28) 

where {pu} is the four-event uniform probability distribution, {pr} denotes the 
thermal probability distribution associated with all possible values for the total 
anyonic charge: total plaquette flux in R trivial and total electric charge in R 
trivial, etc, and where -D({p"^}| |{p^}) — — J2j P] /p'j) denotes the KuUback- 
Leitner pseudo-distance between two probability distributions {p^} and {p^} 
[?7] . We can see that lim^^^oo 7oo(/3) — 0. We conjecture that the formula ( p8| 
is valid in general, with the probability distribution induced by the quantum 
dimensions: 

Pu{q) = (29) 

This relation has been verified for all Abelian quantum double models. 

The other possibility, when studying the asymptotic behavior of the topo- 
logical mutual information, is to let the size of the torus and the size of the 
region R, grow at the same rate. Let i^L x lyL denote the area of the region 
R {v < 1). In that case, keeping 6^ fixed is the only way to make the limit 
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meaningful. Using the fact that liniL^oo^^ = 1 in that case, one finds again 



Eq.(28l, but the probabilities are slightly different now. For example 



Pl^oniR) ^ ^ 2(1 + 9^') ' ^ p:vc„(^)- (30) 

These quantities are still the probabilities corresponding to the value of the 
total charge (flux) sector for region R but subject to the global flux (charge) 
neutrality condition. 

Eqs. ( [28p0l ) are very interesting in that they allow to extract a scaling law 



for the topological mutual information. In the simultaneous limit, i.e. for a fixed 
value of u, the topological mutual information only depends on the temperature 
and size through the parameter t = tanh(/3)^ . In particular, a fixed value of 
t, and thus a fixed value of the topological mutual information, corresponds to 
the following relation between size and temperature 

Pit, L) = In L - ^ ln(i In ^) + 0{L-^), 



dT{t,L) _ -1 

dL L{\nL-\H\\n\) + 0{L-^)Y^^ ^ ' 

For a general Abelian quantum double, based on a group G — I^Ki x • • ■ x^it^, 
these formulae generalise to 

Pit, L) = InL - i ln(-i In ^) + 0{L-'), 

where K = Ki . . . Kj., and where the scaling variable is now defined as 

1 - e-'3^^ . L= 



t - 



(r 



+ {K -l)e-fi^El ' 

where Ai? denotes again the energy associated with an excitation. 

These relations tell us how an increase of the size of the system should be 
compensated by a decrease of temperature in order to maintain a fixed value of 
the topological mutual information. We believe that they constitute a qualita- 
tive nuance from the fact that at any fixed finite temperature, the topological 
mutual information, or the topological entropy, asymptotically vanishes when 
the size of the system is increased [15 . In particular, they show that the rate at 
which the temperature should be decreased, in order to maintain a fixed value 
of 7, decreases with the size of the system. 



3.4 Mean value of the ^-matrix 

As we have seen, the sub-leading term in the area law for Ir-.rp contains entropic 
information about the thermal probability distribution of topological sectors in- 
side region R. In fact, we now wish to argue that this probability distribution is 
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the fundamental quantity controlling in principle (up to implementation-specific 
problems) the usefulness of topological quantum memory and quantum comput- 
ing by anyon braiding. 

For concreteness we take the toric code model that has been initially pro- 
posed as topological memory [3 . Consider the process of creating a pair of 
electric and a pair of magnetic defects, braiding between one particle of each 
pair and subsequently annihilating both pairs, as shown in Fig. [5] The crucial 
property of these anyons is that the operation C/Hopf effecting this process on 
any state of the ground level multiplies it by —1; indeed, t/nopf = ~Ufi^U^^, 
where C/^ and C/jf are the operators describing the processes associated with 
the electric and magnetic anyons separately, i.e., unlinked. Thus, the thermal 



Figure 5: Braiding of anyons in the toric code. A pair of electric defects 
(spheres) and a pair of magnetic defects (cubes) are created. One electric and 
one magnetic defect are wound around each other, then both pairs are anni- 
hilated. The anyon trajectories define a two-dimensional projection of a Hopf 
link, whose components enclose regions Ri and i?2. 

expectation values are immediately related to (marginals of) the charge proba- 
bility distributions since [/^^ and U^^^ measure the electric and magnetic charges 
inside their regions: 



(C/Hopf)/3 - - {U^RjpiU'^'Jp = - (p:ve„(^l)-Kdd(i?l))(P?vc„(^2)-pLd(fi2)) • 



As a consequence, the expectation value of i/nopf is controlled by scaling 
variables and expect this will be the case for many anyonic 

models, where similar expectation values control the visibility of interferometry 
experiments [25] . Note that Hopf-link-like processes define the elements of the 
topological ^-matrix and twisted self-braiding of anyons yields their topological 
spin. Thus, these fundamental quantities of the anyon model are degraded at 
finite temperature at a rate controlled by the thermal charge probability distri- 
butions. The latter becomes the object that determines the appropriateness of 
the system to perform quantum computation at finite temperature. 

In particular, /*°p° measures the Shannon entropy of this distribution. This 
is why we believe that it is a good topological order parameter. Note that in the 
ground level, the probability of finding anyons in a given region vanishes, there- 
fore the probability of a given sector q, p'^{R), becomes Sq^i and the distribution 








(32) 
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has zero Shannon entropy. For high temperatures the distribution approaches 
(29), which is dictated only by the quantum dimensions of the anyons. Our 
conjecture (28 1 implies that this is the maximal Shannon entropy available to 
the distribution. 



3.5 Topological Entropy 

We now analyze the behavior of 7' introduced in Eq. ([3| . We start by observing 
that if one writes 7' as Y.R (^rSr, where R G {A, B, C, AB, AC, BC, ABC} (see 
Fig. [1]), then 

Y,^RiNpiR)+N.4R)) = l. (33) 



With this lattice relation and Eq.(22l, one can see that in the limit where the 
sizes of all regions diverge, 

7'->ln2--^-ln(l + e"2/5). (34) 

At fixed temperature, 7' does not vanish as the size of the system grows, a 
behavior that contrasts with that of /topo- It actually becomes independent of 
the system size. Since, at finite temperature, the von Neumann entropy of a 
region is no longer a measure of its correlations with the rest of the system, it 
is not clear whether 7' actually still probes topological order. The discrepancy 
between the behaviour of 7 and that of 7' is made obvious here because magnetic 
and electric defects have been treated on an equal footing right from the start 
{J — J' in Eq.([7])). Our results therefore do not contradict those of [IS]). Note 
that 7' vanishes in the limit where P tends to zero, as expected. 



3.6 The Case of One Plaquette 

In this subsection we develop a quantum simulation of a minimal toric code 
at finite temperature [25!. We therefore consider a four qubit GHZ state, cou- 
pled to a single ancillary qubit, that plays the role of the environment. For a 
particular time-dependent coupling between the system and the environment, 
it is possible to reproduce the exact behavior of topological entanglement as 
a function of temperature, where the latter is now represented by time. This 
can be also viewed as a purification protocol of the topological thermal states. 
An experimental verification of this topological behavior could be demonstrated 
with state-of-the-art technology. 

Consider a single plaquette of the toric code model [23]. The corresponding 
Hamiltonian can be given by 

where we have omitted the term a\a\ as it is superfluous for generating a 
topologically ordered ground state. The corresponding ground state (vacuum) 



14 



is the GHZ state 



le) - \GHZi) 



V2 



(10000) + 11111)). 



(36) 



When the plaquette is occupied by an electric charge, the state of the system is 
|e) = CTilO = ^d'^'-''^*-') ~ 1 11 11))- It is possible to consider the thermal state 
of the system at temperature T = 1/(3, 



-HP 



Xr(e-ff/3) ■ 



(37) 



We can calculate the topological mutual information for this system. Consider 
the splitting of the four qubits in the regions A, B, C and D, as seen in Fig. [6] 
If we plot /topo as a function of the temperature for this system, we get Fig|7] 



i 

A I 



2 \ 



D 



.\ • 

v4 I 
% ^ _ _ 

^'i ■ 

'3! 



I 




Figure 6: (Left) The one plaquette system of the toric code, which is in the 
GHZ state. Each subsystem A, B, C and D consists in a single qubit. (Right) 



The same system without the Hamiltonian ( 35 1 coupled to the environment 



here taken to be one extra ancillary qubit. The latter can be coupled to all the 
plaquette qubits or equivalently to just one of them. 



Since all the interaction terms in the Hamiltonian ( 35 1 commute with each 
other it is easy to evaluate the density matrix. For simphcity we focus on the 
case where J ^ 1 which reduces the density matrix to 

= ^(1 + CTj^CTfafaJtanh/3). (38) 

We would like to simulate the same behavior, but without the background 
Hamiltonian ( |35| . We can reproduce the density matrix, pp, with four non- 
interacting qubits initially in state |^) coupled to an environment. We take the 
environment to be an ancilla initially in state I-^q) = (|0) + |l))/\/2 coupled to 
the system through the interaction 

Hint - ^(T^a:. (39) 
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Figure 7: Topological mutual information, of the plaquette system as a function 



of the temperature, T. When the Hamiltonian (35 1 is present we observe a 
resilience of the topological character (7 ~ 2) for small temperatures, which is 
lost for higher temperatures (7^0). 



Note that an interaction Hamiltonian that couples the ancilla symmetrically 
with all the qubits would give the same results. For U{t) — exp{—iHintt) the 



time evolution due to ( 39 1 gives a reduced density matrix for the four qubits of 
the form 

Pt^tra{Um,i^a){^,i^a\UHt)). (40) 

For particular choices of a time dependent ui, it is possible to make pt identical 
to pf3. Now, the time evolution of pt is identified with the increase in the 
temperature, T, of p/j. Indeed, if we choose the coupling as 

u^{t) = ^ cosh-' J, (41) 



then the temporal evolution of the system without the Hamiltonian ( 35 1 is 
equivalent to constantly increasing the temperature of the same system in the 
presence of the Hamiltonian. 



Thus, the toric code plaquette with Hamiltonian (351 in a thermal state at 
finite temperature, T, has the same topological mutual information as a plaque- 
tte initially prepared in |^) and coupled to an ancilla in state 1+), where time 
plays the role of temperature. By employing larger states and more ancillae, 
one could reproduce the finite temperature topological behavior of larger toric 
code systems. 



4 Non-Abelian models 

We now turn to a class of models featuring non-Abelian anyonic statistics [3]. 
We call these models "non-Abelian superconductors" , as in JJi\ because of the 
braiding properties of their excitations. We mainly focus on a particular model 
based on the quantum double 0(83), because it is the simplest in the family 
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we are considering. This model is paradigmatic though; our analysis can be 
straightforwardly generalized to models defined through the quantum double of 
any finite group. The 0(83) model is also important from a quantum informa- 
tion perspective because it allows to perform universal quantum computation 
[TU] . We start by briefly reviewing some properties of the hamiltonians we are 
considering. Then we show how the entropy of a region can be calculated. This 
result is employed to study the behavior of the topological mutual information 
at finite temperature. 

4.1 The 0(83) superconductor and its eight sectors 



Figure 8: Tiled sphere. Equal numbers refer to identified sets of edges. 

Consider an oriented lattice, A, on a sphere such as the one represented 
in Fig. [8] Let G denote a finite group, where e refers to its neutral element. 
With each edge of A, we associate a |G|-dimensional Hilbert space, !K, with an 
orthonormal 'computational' basis labelled by the group elements {\g) : g G G}. 
We now borrow a series of definitions introduced in We introduce the 
operators L± and T± , which act on the computational basis as L+ (17) \g') = \gg'), 
T+{g)\g') = 5gg,\g'), L.{g)\g') = \g'g-^), T.{g)\g') = 8g-.,,\g'). We also define 
-^^(j) = ^-{9) if s is the origin of an edge j, whereas L^{j, s) — L+{g) if s is 
the endpoint of j. Also, if a plaquette p is at the left (resp. right) of an edge j, 
we define Tg{j,p) — T_{g) (resp. Tg{j,p) — T+{g)). We further define 



E 



Tg,S^m,p). (42) 



je*{s) 



where in the definition of Bg{s,p), ji, . ■ ■ ,jm are the boundary edges of p listed 
in counterclockwise order, starting from and ending at some vertex s. The 
operators Ag(s) and Bg(s,p) commute when they share no link. Otherwise, 
they satisfy the following relations [3]: 



A A 



Bgi Bg^ 



^9192 Bgi 1 



^91 B92 



^919291^^31^ 



These commutation relations represent those of an algebra called the quan- 
tum double or Drinfeld algebra J29j that we denote as D{G). 
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From the operators 

]7^iJ2Ms), and Bp^Beip), (43) 
' ' see 

it is possible to construct a Hamiltonian made of local mutually commuting 
terms, which has the same form as for the toric code, Eq.([7]). 

An elementary excitation of H lives on a vertex and an adjacent plaque- 
tte that we call site. Each excitation is associated with a pair of data: an 
equivalence class of G and a representation of the normalizer of an arbitrary 
element of this equivalence class [23 El US] . Each such pair of data corresponds 
to an irreducible representation of D(G). Quasiparticles associated with the 
trivial equivalence class {e} are called pure charges. Quasiparticles associated 
with trivial representations are called pure fluxes. The rest of the quasiparti- 
cles are called dyons. For the -0(5*3) niodel, the quasiparticles come in eight 
varieties that we label A, . . . ,H [T^]. One can construct operators that project 
the links related to a site onto a deflnite quasiparticle state. This is done upon 
observing that the operators Bh.{s,p)Ag{s) form a reducible representation of 
the quantum double D^S^). Each definite quasiparticle state corresponds to an 
irreducible representation of D^S^) contained in this reducible representation. 
Therefore, the projectors onto a given quasiparticle type can be constructed 
using characters of the corresponding irreducible representation [30l I31j : 

n= E X<}{h,g)Bhis,p)Agis), (44) 

where the characters Xq associated with an irreducible representation q can be 
calculated |ST] . 



4.2 Pinning Quasi-Particles 

We have not been able to fully characterise the spectrum of the Hamiltonian ([t]), 
with As and Bp defined as in Eq.(43). However, we can argue that we actually 
do not need to as far as we are interested in the topological mutual information. 
Excited states of H are tagged by vertices s on the lattice that violate the 
condition (As) = 1 and plaquettes p that violate the condition (Bp) = 1. We 
restrict to that part of the spectrum of H such that excitations are pinned 
at fixed, non-adjacent locations (see Fig. [9] for a representation of a planar 
region). We expect that the more the noise undergone by the system, the less 
the topological mutual information. Considering only that restricted part of the 
spectrum of H can be understood as additional error correction, where some 
plaquettes and vertices are over-protected so that they never get excited (or 
only with vanishing probability) . It is actually easy to write a local Hamiltonian 
that automatically does that. Therefore, we believe that the topological mutual 
information of the modified model can only be larger than that of the full 
spectrum. 
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Figure 9: Square lattice associated with a quantum double model. Extra edges 
closing the lattice onto a sphere are not represented. Thick strokes refer to sites 
of the truncated model where excitations, living on a vertex and a plaquette 
(site), lie apart from each other. 

We can also provide a renormalisation argument to support the idea of work- 
ing with a truncated spectrum. Fixing the distance between two neighbour exci- 
tation sites to, say, two lattice spacings, and increasing the size of the lattice, i.e. 
the number of links, we see that the truncated model has a continuum limit, 
where excitations are pointwise anyonic particles (or superpositions thereof). 
We therefore expect that the potential differences between a situation where we 
consider the full spectrum and a situation where we consider only the truncated 
spectrum vanish in the limit of large lattices. 

Our third argument is that in the case of the TLi toric code, it does not make 
any difference whether we work with the full spectrum or with the truncated 
spectrum. This is most easily seen by considering the partition function of 
this model. Had we computed it for a system with more edges but still only 
I? possible positions for the charge (resp. flux) defects, we would have found 
exactly the same expression. Similarly, one easily checks that the topological 
mutual information remains unchanged. 

4.3 The fusion space and partition function 

The space of n excitations pinned at fixed sites has the structure [3] 

(45) 

where each index qi runs over all possible quasiparticle types A, . . . , H. This 
structure merely reflects the fact that different excitation patterns lead to or- 
thogonal states. Each space ^q^...q^ further splits as 

^q^...q„ ^Xq^® ...(g)Xq^®Mq^...q,^. (46) 

The factor spaces 3Cq. have a familiar meaning. They correspond to the local 
degrees of freedom of the excitations. The point of topological models is that the 
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extra piece A^g^. ..<;„, the fusion space, may have a non-trivial dimension. The 
particles can be fused. That is, two particles can be converted into one using 
a ribbon operator to connect them [3]. In the simple Abelian case discussed 
previously, the fusion rules were simple. For example, two plaquette excitations 
can only be fused in a way that makes them disappear. In a general anyonic 
model, the fusion rules read 

qi, = J2Kb'lc, (47) 

c 

where iV^j, are non-negative integers. We assume that physical states have to 
fulfill some neutrality conditions, i.e. they should be such that fusing all the 
particles yields with certainty the trivial particle, denoted as 1. This assumption 
is legitimate since we work with a sphere [7] . An important result of the theory 
of anyons is a formula for the dimension of A4q^...q^ [T3]: 

dim Al,,...,„ - E • • ■ E ■ ■ ■ <-.9„- (48) 

bi b„-2 

This formula has a structure that makes it easy to compute for arbitrary n. 
dim Mqj^...q^ is the contraction of a (quasi) translationally invariant matrix 
product state f37. The dimension of 3-C[n] can also be calculated easily. Let dq 
denote the dimension associated with the quasiparticle q. For D^Ss) anyons, we 
have that (dyi, rfs, rfc, f^D, rf_E, c^_F, ^G, di/) = (1, 1, 2, 3, 3, 2, 2) [13]. Defining the 
elements of a matrix M as — J2-y N^^dj,dim'K[n] can be expressed as 

dimJ{[7i]= ^ dimJCg,...g„ ^^d^M^ji. (49) 

qi...qn a 

The partition function of the model on a sphere, with n well separated sites 
reads 

ZiP, n)^J2---Il e-^^('-^")d,, . . . d„.<V<,3 • ■ • K^^i. ' (50) 



where E(qi . . . g„) is the energy associated with a configuration qi . . . g„. Just 
like dim!K[n], it can be computed efficiently upon diagonalizing a matrix M, 
whose size is independent of n. 



4.4 von Neumann entropy at finite temperature 

Let us start by considering a scenario where a pair of anyons \qq) is created in 
such a way that anyon q lies in some region A and anyon q lies in the comple- 
mentary region. In this configuration, the von Neumann entropy of region A 
reads [21] 

S^^^"{pA)^S{p^/-) + logdq. (51) 
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The entropy of a region when the system is in a thermal state can be com- 
puted once we are able to calculate the entropy of a region when the system is an 
arbitrary defect configuration. In turn, the latter entropy reduces to computing 
the entropy when there are only two anyons in the system, and one lies inside 
the region we are interested in. 

Consider a tiled sphere as the one indicated in Fig. |8] and let us divide 
it into two simply connected regions, A and B. Let Na and Nb denote the 
number of sites contained in region A and B respectively {Na + Nb = n) . Let 
^A = Qi ■ ■ ■ QNa label the types of anyons living on the sites contained in region 
A and q'^ = q[ . . . q'b label the types of anyons living on the sites contained in 
region B. The total Hilbert space can be decomposed as: 



nn] = 00 ^Cq, ® scq', ® -ML ® ® Ki- 



rn 



qv 



(N.B. The spaces -Ml^ are one-dimensional.) The decomposition (52) induces 
the following representation of the thermal state of the 0(83) superconductor 

Ah =000 0- 

q.A q^ b=A Ml = l M2=l SA = 1 sb = 1 



Z{f3,n) 



\qA b,^ii,SA;q'B b,H2,SB){qA ^ b, ^ii, SA;q'B ^ 6,M2,ss|, (53) 

where E{qA) = Y.j'A^ilj)' d,{qA) = YlfAdiQj), A^i (resp. fi2) denotes the 
possible channels through which the anyons (resp. q'^) can be fused into 
b (resp. b), and sa (resp. s^) is a collective index for the internal degrees of 
freedom of the quasiparticles contained in region A (resp. B). 
Defining 

ZA{P,b) ^ ^d(q^)e-''^('i-)dimA^^^, 
qA 

and ZB{f3,b) likewise, we get that the reduced state of region A reads 



dimA^° 



di<iA) 

p^-0 



e'PEiqA)Zg{(3^b) 



trslqyi -> b,fj.i,SA;b){qA -> 6,^i,sa;6|. 

(54) 



It is worth observing the resemblance between Eq.(54l and the reduced density 
matrix for the toric code. From Eq.(54) and Eq.(20l, we deduce an expression 
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for the entropy of system A: 



^^P^^-- 2. ^ip;;^) i^^i ^p;^) J-^? (P^)\ 

qA,b,)J.i,SA 

(55) 

Just as for the partition function, S{pa) can be computed by diagonalization of 
an auxiliary matrix. Finally, the von Neumann entropy of the whole sphere can 



be efficiently computed from the partition function using again Eq.(15). 

These expressions are valid for any quantum double model on a lattice such 
as the one depicted in Fig. [9] Remarkably, the only microscopic data involved 
in the expressions for S{pa) and S{pth) are the tensor N governing the fusion 
rules, the energies associated with each quasiparticle type and their internal 
dimension. It is therefore tempting to conjecture that these expressions are 
valid for other anyonic theories. 

The fact that the fusion rules tensor N*.^, appears in the expression for the 
finite temperature entropy of a region marks a qualitative difference with re- 
spect to the zero temperature (ground state) case. This is interesting in that it 
shows that entropic quantities could be used in order to discriminate between 
different topological models. For example, at finite temperature, the quantum 
double models based on Zg and S3 exhibit different region entropies, while these 
quantities coincide at zero temperature. We understand this difference as fol- 
lows: at finite temperature, quasiparticle pairs can be spontaneously created 
from the vacuum and can move around on the lattice. When they braid, they 
give rise to qualitatively different statistics, which in turn give rise to different 
values for entropic quantities. 

4.5 Numerical Results 

The results derived in the previous section, and a rigorous calculation of the 
ground state entropy, shown in Appendix [Clhave allowed us to study numerically 
how the topological mutual information behaves as a function of /3. The systems 
we have considered are four tiled spheres, all with 96 x 96 plaquettes. The first 
sphere contains 64 sites, the second 144, the third 256 and the fourth 576. 
Although these systems are small when compared to their Abelian counterpart 
we have considered on Fig|4] they are large enough to show that non-Abelian 
systems are affected by temperature in the same way as the toric code [12 . 



Fig. 10 is presented to analyze the existence of scaling laws in the non- 
Abelian case. There, the topological mutual information is represented as a 
function of nexp(— 2/3). We can see that the curves tend to collapse when the 
number of sites is increased, that is, the distance between the four curves re- 
duces. This fact supports the idea that for large systems, the topological mutual 
information only depends on the 'volume' of the system and the temperature 
through the product nexp(— 2/?). This result is consistent with the scaling re- 
lation found for the toric code. 
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Figure 10: Topological mutual information as a function of nexp(— 2/3) for the 
DiSs) model. 



5 Conclusions 

In conclusion, we have studied how topological phases of lattice systems behave 
as a function of temperature, using the topological mutual information as an 
order parameter. We have focused on two specific models; the toric code and 
the 0(83) quantum double model. These systems turned out to be afi^ected by 
temperature in much the same way. We have studied how the topological mutual 
information depends on size and temperature through a simple scaling variable. 
In particular, we have shown how an increase of system size can be compensated 
by a vanishing decrease of temperature. It would be interesting to see whether 
our findings are also valid for other systems, such as a FQHE sample for instance. 
We have also seen how temperature could allow to resolve between models whose 
topological entropy are equal at zero temperature. We have also introduced 
a class of simulations where running time simulates increasing temperature. 
Though exemplified by a minimal plaquette in the Abelian model, it can be 
extended to more complicated situations. We have also discussed why, at finite 
temperature, hopo should be used as an appropriate topological order parameter 
instead of the topological entropy. We have shown how /topo relates directly to 
the thermal probability distribution of topological charges, and how in turn 
this distribution is related to the mean value of ^-matrix elements and controls 
the visibility of anyonic processes, and therefore the usefulness of topologically 
ordered systems as quantum memories. 
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A Some calculations for the toric code 

Most of the calculations presented in this appendix are based on repeated use 
of three elementary identities: 



fe=0 ^ 



k^O ^ ^ k^O 
k even 



y(^)ka''^a^il + a)^^N a(l + a)""-'. 
t^,\kj da 

One finds that the partition function reads: 

= [(2cosh/3)^' + (2sinh/3)^']2, (56) 

as derived in (TSl [H] . 

Next we want to compute the quantity 

W{P, R,L) = -Y, 4C(qK) ln(4C(qfl)), (57) 
qi? 



appearing in (21 1. For that, it is convenient to consider plaquette and vertex 
excitations separately: = (tpji^cji). Let Np{R) (resp. N^,{R)) denote the 
number of plaquette (resp. vertex) excitations whose edges are all contained in 
R, and let us define N^{R) = - Np{R) and Nl{R) = L^-N^{R). With these 
notations, C(q_R) can be decomposed as 

C{cIr) - J'(|0«|,iV;(i?),e-2'^) n\cR\,K{R),e-'^), 



where 

N 



7{m,N,a)^ ^ (^)a- = Uil + af + {-ir{l - a^}. (58) 



n-\-m even 
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In order to get an anlytical expression for W{l3, R, L), it is sufficient to evaluate 
sums like: 

S,{R) = ^e-2''l*«lT(|0fl|,7V;(i?),e-2^), 
<i>R 

and 

S2{R) - e-^^^^^^y{\<l>RlN;,{R),e-'f') In 5'(|</>fl|, iV;(i?), e^^/j)^ 
S,{R) = ^e-2;3|0«|j>(|^^|,Ar^(i?),e-2/3) 

^ m=0 V / 

= + (59) 

Splitting 5*2 (i?) as a sum over even values of and a sum over odd values of 
\4>r\, we get that 

J{(l + e-2/3)A'p(«) + (l_e-2/3)A^p(fl)}{(l + e-2/3)Ar;(fl)^(^_g»2/3)iV;(fl)| 



ln[i{(l + e-2/3)<(«) + (1 - e-2^)^^(^)}] 



'2 



+ + e-2/3)^p(«) - (1 - e-2'3)^p(«)}{(l + e-2/3)^;(«) _ (i _ e-2/3)^;(«)} 



4 



ln[i{(l + e-2/3)<(«) - (1 - e-2/')^;(«)}] (60) 



With these results, we find Eqs.( p2p3 l 



B Rough edges and smooth edges 

When dividing a lattice into four regions as indicated in Fig. [l] one has to be 
careful that the links on the boundary between two regions should be attributed 
only to one of them. On FigjlT] we show how we have chosen to divide Fig{T] 
Data related to the division are collected in Table 1. 

In the Abelian case, the entropy of a region R in a zero temperature thermal 
state (i.e. in an equal mixture of four basis states of the ground states) is given 
by (|9_r| — 1) In 2 — e^. The cr correction appears because we have worked with 
a model defined on a torus, and takes into account the fact that the region R 
might be contractible. The "area" \dR\ is given by the total number of crosses 
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minus the number of crosses fully inside R minus the number of crosses fully 
outside R. 

The non-Abelian 0(83) case has been studied on a sphere. There the ground 
state entropy is given by — l)ln6, see Appendix C. As shown in Fig. [sj 

there is a region within the sphere made of square plaquettes only. (We have 
chosen our tiling of the sphere in this way.) Let L denote the (linear) size of this 
region, and k denotes the linear size of region A. (All triangular plaquettes are 
assumed to be inside region D.) In order to ease the computations, we choose 
L to be a multiple of fc, and we impose that the density of sites, s^, is such that 
is a divisor of k. (Triangular plaquettes are assumed to support no site.) 





Figure 11: Division of the lattice into four regions used in numerical calculations. 
Delimitations of a region by dotted lines indicate that the corresponding links 
do not belong to the region. 



Region 


"Area"-1 


e 






"Volume" 


A 


4fc- 1 


-log 4 




{k-lf 


fc2 


B 


4fc - 5 


-log 4 


P -Ak + A 


{k-lf 


fc2 


C 


6k- 1 


-log 4 


2fc2 - k 


(fc-l)(2/c-l) 
- 4fc2 - Ak 


2fc2 


D 


8fc- 2 





- Ak^ - 6k 


L2 - 4fc2 


AB 


6fc- 3 


-log 4 


2fc2 -3fc + 2 


(fc-l)(2/c-l) 


2fc2 


AC 


8fc- 2 


-log 4 


3fc2 


-Ak + 2 


3fc2 


AD 


8fc- 3 





- 3k^ - Ak 


-3k^-Ak + l 


L2 - 3fc2 


BC 


8fc- 3 


-log 4 


3fc2 - 4fc + 2 


'dP - 4A: + 1 


3fc2 


BD 


8fc-2 





L2 -3fc2 -8fc+ 1 


£2 - 3A:2 - 4/c - 1 


L2 - 3fc2 


CD 


6fc-3 





- 2P - ik 


£2 - 2A:2 - 3/j + 1 


L2 - 2fc2 


ABC 


8fc-2 


-log 4 


4fc2 - 2fc + 1 


{2k - 1)2 


4fc2 


ABD 


6k- 1 





L2 - 2P - 5k 


L2 - 2A:2 - 3A; - 1 


L2 - 2fc2 


ACD 


4fc-5 







L2 - fc2 _ 2fc + 3 


L2_fc2 


BCD 


4k- 1 





-k^ - Ak 


-k^ - 2k - 1 


L2_A:2 



Table 1: Some data useful to compute the von Neumann entropies of various 
regions of a tiled torus at finite temperature. The regions are those depicted in 
Fig. 11 e labels the correction to the ground state entropy due to the topology 



of the considered region. Np (resp. N^) denotes the number of plaquettes (resp. 
crosses) fully contained in the region. 
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C Ground State Entanglement of the D{S^) su- 
perconductor 

Here, we generalize some of the calculations presented in [33| and derive a simple 
exact expression for the entropy of the reduced matrix pA of the ground state 
If) in a region A. {B here denotes the complementary of A). We need to 
distinguish amongst three types of vertices, those in A, those in B and those 
touched by edges belonging to A and edges belonging to B. The set of vertices 
of the latter kind will be referred to as dA. We consider only regions with the 
following properties 

1. Both A and B are connected in the sense of vertices. That is, every two 
vertices inside A (resp. B) can be connected by a path of adjacent vertices 
inside A (resp. B). 

2. For each s e dA, there exist s' £ A, s" € B adjacent to s. 

The (unnormalized) ground state of the D{S^) model defined on a sphere 
can be written as 

10= ^9i(«i)---^s«(«Jv)|e---e) = ^ 7r(g)|e---e), 

where g = (gi, . . . ,5Ar), G = ^3 x • • • x 5*3 and tt is the representation of G 
defined by 7r(g) = Ag^{si) ■ ■ ■ Ag^{sN)- The action of 7r(g) can be decomposed 
as a product of unitaries acting on the edges of the lattice: 

7i"(g) = U{s,.s,){g^,gj), 

where the operator J7(s.,sj) acts as follows on the 'computational' basis states 
{\x) : X e 53} of the edge {si,Sj): Ut^si. sj)\x) — iQi^gJ^) if the edge goes from 
the vertex to Sj. Therefore U(^s. s ){giTgj) acts trivially on an edge {si,Sj) 
only if gi equals gj and belongs to the centraliser of S3, which is {e}. Therefore, 
the representation tt is faithful and we can (and will) identify G with 7r(G). 

We now borrow some definitions from |33j . We define two normal subgroups 
in G: 

Ga = {g e G\gj = e if ^ A}, 

GB = {geG\g,=eiis,^B}, 

and we call Gab the quotient group G/{Ga x G^). 

Any element of Gx {X — A, B) acts only on X and we have the following 
partition of G: 

G = U[h]eGAB{(gA«igi3)h|gA e GA,g6 e Gb}, 
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where we can fix any representative h for the class [h]. This decomposition 
impHes that the ground state can be expressed as 

\0 ^QA®QBY.h\e---e), (61) 

[h] 

where Qx = Eg^sGx S^' 

For our purposes, we need to quotient further Gab- So, we consider the 
diagonal of G, Gd — {r — {r, - ■ ■ ,r)\r £ S3}, and the quotient set, G^s/Gd, 
defined by the (right) equivalence relation: [h] ~ [h'] iff 3 r € 6*3 s.t. [h] ^ [^'r]. 
Clearly, jGAs/Gdl = {GabI/ISsI We denote [[h]] the elements of G^s/Gd- 
Since Gd acts trivially on |e . . . e), |^) can be further decomposed as 



\O^Qa(S>Qb Y1 ^^r\e...e) = l^alQA ® Qi3 ^ hs|e . . . e), (62) 
reS3 [[h]] [[h]] 

with h = h.A'E)h.B, and where (resp. h^) does not necessarily belong to Ga 
(resp. Gg). In order to conclude our calculation of S{pa), we use the following 
lemma. 

Lemma 1 Let g denote an arbitrary element of G and let g = gA ® gs denote 
its decomposition into an operator acting on A and an operator acting on B. If 
for all r e S3, [g] 7^ [f], then (e . . . e|gx|e ... e) = /or both X = A,B. 

Proof: Let us suppose that it is false, [g] ^ [f] iff there are at least two 
vertices Si, Sj e OA such that gi =/= g.j. Let us consider a path Tx connecting Si 
to Sj through edges within X. Since 

(e...e|gx|e...e) = J]^ {e\g-^gp\e), 

we have that (e . . . e|gjs:|e . . . e) 7^ only if ga — gp for all pairs of adjacent 
vertices Sq, sp. Thus (e . . . ejgjfle . . . e) 7^ only if gi = g^ for all vertex Sa G 
Fx, and in particular gi = gj, which is the desired contradiction. 

This lemma allows us to prove that the states {(h^ 'S) h.B)\e . . . e) : [[h]] €E 



G^s/Gd} form a bi-orthogonal set, so that the expression (62) is actually 
a Schmidt decomposition of the ground state. If [[h]] and [[h']] are different 
equivalence classes, then (e . . . ejh^^^hj^le . . . e) = for both X — A,B. Indeed, 



[M ^ 



of Ref. 



[h']] means that [h] ^ [h'r] for all r e S3. That is [h-^h'] [f] and the 



lemma T] shows that the desired property holds. At this point, the construction 



can be used to first prove that the set {{QA^i-A ® QBii-B)\e . • . e) 



[[h]] <E GAB/Gd} is also a bi-orthogonal set, and then on to show that 

S{pa) = log2 = log2 \S3\{NaA ~ 1), (63) 
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where Nqa is the number of vertices on the boundary of A. 

Remark: The argument just presented for 6*3 is actually valid for any finite 
group G. If G has a non-trivial centre Z{G) ^ {e}, then one would merely start 
the construction from the group G/Z(G)d instead of G. 
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